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Abstract 

A new mechanism is proposed to explain the "anomaly" in Z — > bb resulting 

in a prediction of a new, non- sequential fourth family whose masses could all 

be below Myy, thus opening up an exciting prospect for near- future discoveries 

at LEP2 and possibly at the Tevatron. 
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Precision tests of the Standard Model (SM) have reached a level where it "might" now 
be possible to look for indirect evidence of new physics and/or new degrees of freedom. One 
example is the apparent discrepancy between theory and experiment in the value of the ratio 
Rb = T(Z — ► bb)/T(Z — ► had). This discrepancy which increases with mt, reaches the 2 a 
level when m t reaches 175 GeV. If one also includes the apparent disagreement between the 
QCD coupling as determined at "low" energy and evolved to Mz with that determined by 
the Z-lineshape, one is tempted to think that one might be already seeing some new kind of 
physics. It is therefore very crucial to confirm or disprove these so-called discrepancies. Let 
us nevertheless assume that they are there to stay and examine what kind of new physics 
that can be possible and what predictions that can be tested in the near future. (Even if 
the discrepancy were to disappear, this would put a severe constraint on this type of new 
physics.) 

In this letter, a mechanism is proposed to explain the apparent increase of Rb and to make 
further predictions on other branching ratios, and ultimately on the new physics concerning 
the mechanism itself. This mechanism is based on the assumption that there is a new, heavy, 
non- sequential down quark (Q = —1/3) ( part of a new family) with mass greater than 46 
GeV and whose qq bound state(s) mixes with the Z boson. By non-sequential, we mean that 
the fermions of the new family does not have ( or has little) mass mixing with fermions of 
the other three generations. Other than this being a working assumption, a realization of 
this scenario is given at the end of the paper. 

For this paper, we shall quote a few relevant observables |1],0: F z (GeV) = 2.4974 ± 
0.0027±0.0027(2.496±0.001±0.003±[0.003]),i4 = T(bb)/T(had) = 0.2202±0.0020(0.2155± 
± 0.0004), R c = Y(cc)/Y(had) = 0.1583 ± 0.0098(0.171 ± ± 0), where the numbers in 
parentheses are the standard model expectations, and R e = T(had)/T(ee) = 20.850 ± 0.067, 
R^ = T(had) /T(fip) = 20.824 ± 0.059. 

Let us denote this non-sequential family by (TZ, V) for the quarks and by (A/ - , £) for the 
leptons. For reasons to be given below let us assume that the (Q = —1/3) quark has a mass 
m-p < m-n- We also assume that the up- type quark 1Z is heavy enough so that 7ZTZ bound 
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states are well above the VV open threshold. The VV bound states can be described by 
Richardson's potential. Such an analysis has been carried out long ago by |3[] for the 3 Si 
ti bound states, but unfortunately in the now-obsolete range of mt ~ 40 — 50 GeV. This 
analysis can however be applied to any quark in a similar mass range or higher, especially 
for our case where m-p > 46 GeV. (The mass shift of the Z boson due to this mixing is 
negligible ||.) 

VV bound states which can mix with Z are either vector, axial vector, or both. In what 
follows we shall neglect the mixing of Z with the axial vector states since it goes like (3 3 j|,[3| . 
Consequently we shall focus only on the vector meson ( 3 Si) bound states. In particular, we 
shall first examine the mixing of the ground state IS" with Z. In the mass range considered 
here, the ground state IS* is sufficiently far from open-P threshold so that the mass-mixing 
formalism can be applied. Denoting the IS" (J PC = 1 ) state by V°, the result of V° and 
Z° mixing is given in terms of the mass eigenstates [[| 

9 9 

\V) = cos-\V ) - sin-\Z ), (la) 
9 9 

\Z) = sin-\V ) + cos-\Z ), (lb) 

for the mass eigenvectors and where 

9 = sin" 1 {5m 2 /A 2 ), (2) 
with A 2 = ^Kr ir vo M v -^q o +ir Zo M Zo r + ^ m 2 )2 ]i/2_ 6m 2 ig the off . diagonal element of the 
mass mixing matrix and is given by || 
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where F v = 2v / 3|*(0)|^/M Vb . The term inside the square brackets represents the vector 
coupling of the V quark to the Z boson. 

Let us assume that My > Mz and since present experiments are carried out on the 
Z resonance, we need only to look at Eq. (|I~b| ) to see how the presence of Vq modifies the 



coupling of Z to "light" quarks and leptons. This, as we claim in this manuscript, is a 
possible source for the discrepancy seen in T(bb). From Eq. flTED, one finds the physical Z 
couplings to a given fermion / to be 

v.A ■ 6 v,a . & V,A /a 
q' t f = stn-q, /f f + cos—q' f f, (4) 

where V and A stand for vector and axial-vector couplings respectively. 

The most obvious source of the coupling of V to // is via 7 and Z and evaluated at 
s = Mf. The electro weak source alone however gives only a small change to Rb and in the 
wrong direction, and this worsens when Vq is close in mass to Z. A new and unconventional 
coupling of V to b ( and to other normal fermions as well) is needed, not only to compensate 
for this small electroweak change but also to bring Rj, closer to its experimental value. To 
this end, let us write 

9l A {f = F v G v As = Ml)+g v n >* (5) 



where 



c v (M 2 ) _ i QfQv 9 2 9j_9v , s 

G f {Ml) = cos 2 6 w M z -m{^M Zo T z ' (6b) 
and where g J' A is the vector (axial- vector) coupling of the Z boson to the fermion /, and 
g^p = —(1 — (4/3)sin 2 9w)/^ and g£ = 1/4. Qf and Qv{= — 1/3) are the electric charges. 
9newj 1S the coupling of Vq to a fermion / coming from some new physics. Since 
and M Zo T Zo ^ 3xl0^ 2 M|, we can safely neglect the 7 contribution in Eq. (|6a|) (it contributes 
negligibly to the present discussion). 

For the mass range considered below (shown in the Figure), namely m-p ~ A8GeV — 
53GeV, |^(0)| is such that \5m 2 \ < |M£ - iT Vo M Vo - Mf + iT Zo M Zo ) 2 \/2 and conse- 
quently 



with cos | ~ 1. Typically, 8/2^2 — 3 x and the deviation of cos| from unity will be 
of order 10 -4 and can neglected considering the present level of precision. 
The modified couplings of Z to a fermion / are now 

gj = (1 + vlw + vlneJgJi (8a) 



gf = (l + rfi w + rif inew )gf, (8b) 
where W stands for electroweak and the r/'s are complex numbers and are defined by 

r,J; A = sin-F v Gf A (s = M 2 z )/gJ> A , (9a) 

vj',new = sin-g^ e i/gJ' A , (9b) 

where the explicit forms for rjj'^ and nJ' A ew can be obtained by using Eqs. (|5|, |Ba] , prj| ,[7D . For 
simplicity, we shall now assume that g^ ew / = 9nem 7^ and g A ew j = so as to reduce the 
number of parameters and to study the implication of such an assumption. The introduction 
of a single new coupling is what we meant by universality earlier. We shall argue below why 
we expect such a behavior. (The inclusion of g A ew j is quite straightforward.) 

In computing the Z widths using Eqs. ( |3a] , |8"bl ) and the range of mass mentioned earlier, 
one can safely neglect terms proportional to (Re rj) 2 and (Im rj) 2 since they turn out to 
be at least two orders of magnitude smaller than terms proportional to Re rj (assuming 
9new f < !)■ (Considering the present level of precision, their inclusion is irrelevant to the 
present discussion.) With this remark in mind, the decay width for Z — ► // is now given by 

r(z^ff) = r s f M (i + 5U, (io) 

where / = q,l and where 



r f _ ^{{g- f njte rjw + ne r, new ) + {gyrne 7jfr) 

new i V\2 j_ ( „A\2 ' v ' 



2(( gj 7(fc i + flei ? L) + {9 A fRe^) 

(gj) 2 + (gfy 

In Eq. (0), Tf 1 contains various radiative correction factors as well as mass factors such 
as defined in Ref. ( @). We find 



T(had) = T^(had) + 5 u new (T^ + I*» ) 

+ s d new (r s d M + v s s M + r s b M ), (12a) 

R = gro+jy fl2W 

where Rf = T(Z —>■ qf(jf) /Y(had). The central theme of this paper is the use of Rb to obtain 
information on the model proposed here. By using Eq. ( |12b| ) for Rf,, one can extract the 
parameter Rer\( new and consequently the common parameter sin^g new as a function of My . 
This will then be used to make predictions on various ratios mentioned above and also on 
the total Z width. In particular, R e fX will be used to constrain the range of allowed My . 

We shall use the following hadron ratios: R b = 0.2202 ± 0.0020, R§ M = 0.2155, R S C M = 
0.1721, R S S M = 0.22, R S U M = 0.1722. The Standard Model predictions |] given here are 
for m t = 170 GeV. A more extensive analysis using the range of m t given by CDF and 
DO will be given elsewhere. Notice that the Standard Model values quoted here are rather 
insensitive to the Higgs boson mass. 

We predict: R c = 0.165 q= 0.003, R s = 0.225 ± 0.002, R u = 0.165 q= 0.003 to be compared 
with R c ,exp = 0.1583 ± 0.0098 (more than la lower than the standard model prediction). 

Notice that an increase in the ratio for a down-type quark corresponds to a decrease in 
the ratio for an up-type quark and vice versa. This happens because Rer]J new is positive for 
f = u,c and negative for f — d,s, b. (V is a VV bound state.) 

Let us turn to R e ^. The experimental values are R e = T(had)/T(ee) = 20.850 ± 0.067, 
R^ = T (had) /T (lip) = 20.824±0.059 to be compared with the Standard Model expectations 
(m t = 170 GeV): R s e M = R™ = 20.774, . . . , 20.754 for m H = 100, . . . , 1000 GeV. Although 
these numbers are consistent within errors (except for R e ) , experimentally there seems to be 
a tendancy for an increase in these ratios. (Even when we take into account the spread in 
m t , -Rfj/ < 20.78.) Our predictions for R e = R^ are shown in the Figure (curves labeled by 
300 and 700) as a function of My and for two different values of m#, namely m# = 300, 700 
GeV. They are obtained by using R^ m i n ( the theoretical curves obtained by using Rb, ma x 
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are out of scale in the Figure presented here). The predicted regions lie above these curves. 
The two vertical lines represent the lower bounds on M V) namely M Vo = 95.6, 103.25 GeV 
for ttih = 300, 700 respectively, coming from the Tz constraint. The arrows indicate that 
the regions allowed by T z are to the right. Finally the two horizontal lines delimit the 
experimentally allowed region which we take to be the overlap between R e and R^. There 
we take 20.883 as the maximum ( from R^) and 20.783 as the minimum (from R e ). We 
only show m H up to 700 GeV to be consistent with the global fit although larger values are 
entirely possible. (For m H = 1000 GeV, the lower bound on My is 95.2 GeV.) 

The allowed regions are the ones bounded by the theoretical curves, the vertical lines 
and R e ,ex P - From the Figure one can see that the allowed region for m H = 700 GeV is 
significantly larger than that for m H = 300 GeV. This implies that if the resonance were to 
be found say at 96 GeV one would infer m H > 700 GeV, while if it were found at 103.6 GeV 
one would have a looser bound, namely m H > 300 GeV. A lower M Vo implies a higher lower 
bound on m H . We can also conclude that, for m H < 1000 GeV, the resonance mass which 
is compatible with all available data is bounded from below by 95.2 GeV. This corresponds 
to m P > 48.5 GeV. 

From R b , we can extract g new as a function of M Vo and use these values to compute T Vo . 
The dominant contribution to Ty () turns out to come mainly from this g new with 7, Z, and 
three gluon processes contributing a small amount to the total width. For M Vo > 95.2 GeV, 
the lower bound on T Vo is found to be 1 GeV and increasing to approximately 8 GeV for 
My PS 104 GeV. The expectation is in general a few GeVs for the width in our scenario while 
a standard heavy onium will have a width of at most a few Me Vs. Notice that the shift in 
width due to the mixing with Z is small with respect to the above intrinsic width. A further 
prediction is the fact that, in our scenario, the new coupling is universal so that Vq couples 
equally to quarks and leptons of both up and down types. This implies that r(V — > 11) = T 
and T(V — > qq) = 3r . The prediction for the branching ratios is Bi = 1/24 and B q = 1/8 
where I — e, /i, r, v e ^ T and q — u,d, s, c, b. 

Let us now turn to the other members of this non- sequential family, the TZ quark and the 
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leptons M and £. This is where the S and T parameters |7| come in. Since this new family is 
non- sequential, there is no reason to expect the mass splitting between up and down members 
to be "similar" to the other three families. We use the results of |I[ (a seven parameter fit) 
which show the allowed regions in the T' new — S' new plane where T' new = T new + Tm h and 
S'new = Snew + Sm h - In view of the discrepancy between the SLD and the LEP asymmetries, 
we shall use as the allowed region the overlap of all data except the SLD asymmetries. In 
particular we would like to ask whether or not all of these new particles can be lighter than 
80 GeV, an exciting scenario since they can all be accessible to LEP2 in that case. Since 
the possibilities are many, a few examples will be illuminating. 

Let us take m-p = 50 GeV and m# = 700 GeV. The Higgs contribution to S and T is given 
by Sm h = 0.045 and T Mh = —0.132. For = 60, 80 GeV, one has S neW)quar u = 0.138, 0.107 
and T newmark = 0.497, 0.686 respectively. For the leptons, we shall assume || that AT is a 
Majorana particle and that its mass ( as well as that of S) is greater than 46 GeV. This is 
the only direct constraint one has on the leptons. We shall use two representative sets of 
values. 1) For = 46 GeV, S new ^ ept on = 0.055,0.027 and T new j epton = —0.004,0.011 for 
m £ = 60, 80 GeV respectively. We then get the following results in terms of {S' new , T' new ) for 
the pair (m n ,m e ). We obtain: (0.238, 0.361) for (60, 60) GeV, (0.207, 0.55) for (80, 60) 
GeV, (0.21, 0.376) for (60, 80) GeV, and (0.179, 0.565) for (80, 80) GeV. 2) For m M = 60 
GeV, S neW) i ept0 n = 0.073, 0.054, T neW)lepton = -0.013, -0.007 for m E = 65, 80 respectively. 
We obtain: (0.256, 0.352) for (60, 65) GeV, (0.225, 0.541) for (80, 65) GeV, (0.237, 0.358) 
for (60, 80) GeV, and (0.206, 0.547) for (80, 80) GeV. All of those values are inside the 
allowed region shown in |J. 

From the (not-exhaustive)examples given above, it is clear that the S and T parameters 
certainly allow for the existence of this new, non-sequential family and that the T parameter 
tends to favor a value of 1Z mass lower than 80 GeV (and in no way should it be more than 
90 GeV). This opens up the possibility that the whole family can be found by LEP2. First 
the R ratio would be 16/3 or at least 14/3 (if the 1Z quark mass is above 80 GeV). Secondly, 
there would be two narrow resonances: the first one being the VV bound state and the 
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second one being the TZTZ bound state. Since these quarks do not mix with with those 
of the other three generations, V will be relatively stable and the search for 1Z will be an 
interesting problem in hadron colliders. Finally there will be an unmistakable signature for 
the charged heavy lepton in LEP2. These issues will be explored elsewhere. 

Finally we would like to say a few words about a possible origin of the coupling g new . 
First solving for g new using Eq. ( |12b|) , one has for example the following range: g n ew,min = 
0.135 — 0.377 for My = 95.5 — 104.5 GeV where Rb,min has been used. (It turns out 
that the constraint coming from R e ^ gives an allowed range for g new ranging from the 
previous g n ew,min to a slightly higher value for each My .) It is clear that the larger My 
is (and consequently less mixing with Z) the larger g new should be in order to preserve the 
"anomaly" in R^. The following discussion will be very speculative but helpful to illustrate 
a few interesting scenarios. 

Let us now imagine there is a four-fermi coupling of the form: (g 2 / 'h 2 )P^ '^P 'j f ', where 
/ denotes any fermion of the first three generations. The translation of this coupling into 
g new is of course highly mo del- dependent. The simplest (and most likely easiest to be ruled 
out) scenario is to assume that the above coupling comes from the exchange of some vector 
boson with point-like coupling to the fermions. We would then identify g new = (g 2 / A 2 )Fy, 
meaning that it can be described by the wave function at the origin, li g 2 /4n = 1 (a 
strong coupling scenario), we get A = 163 — 104 GeV for My = 95.5 — 104.5 GeV, while 
if g 2 /^ = 2.5, we would get A = 259 — 166 GeV for the same range. These scales are 
"uncomfortably" low. Another scenario is to assume that all fermions are composite (for 
instance they could be bound states of a scalar field and a fermion field). For definiteness, let 
us assume that only the fermionic constituents carry color. A four-fermi coupling given above 
would be diagramatically similar to the quark diagram for meson-meson scattering except 
that here we would have a scalar line instead of one of the two quark lines. It follows that 
g new is not necessarily given by the wave function at the origin. We shall assume that we can 
write g new = (g 2 / k 2 )g 2 H Fv where g 2 H represents the rescattering of the scalar components. If 
gl/A-n = g 2 H /A7T = 1, 2.5, A can be found to be respectively 580-371 GeV and 1.45-1.04 TeV 
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for M Vq = 95.5 — 104.5 GeV. Here A would represent the compositeness scale, which as we 
have seen could be in the (low)TeV range. A model which we are currently investigating is 
similar in spirit to the Abbott-Farhi model except that the confining gauge group is not the 
electroweak group but a horizontal gauge group which we take to be SU (2) L ®SU(2) R . In this 
model, there remains a residual global horizontal SU (2) with composite fermions forming a 
triplet (the three standard families) and a singlet (the non- sequential fourth family) under 
that group. A full discussion of the model is beyond the scope of this paper. 

We have presented a simple scenario to explain the "anomaly" in and, as a con- 
sequence, we have made a number of predictions including the presence of a new, non- 
sequential fourth family whose masses could be all below M w , an exciting prospect for 
near-future discoveries. 

This work was supported in part by the U.S. Department of Energy under grant No. 
DE-A505-89ER40518. 
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FIGURES 

FIG. 1. The allowed regions in the R e (= T(Z had)/Y{Z -> ee or fifi)) - M Vo plane. The 
theoretical predictions lie above the curves labeled by 300 and 700 for two different values of 
m,H- The vertical lines with similar labels represent the regions (to the right) allowed by Tz- 
The experimentally allowed regions lie between the two horizontal lines labeled by R e ,exp- The 
intersections of these lines represent the final allowed regions. 
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